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Abstract 

We present a high sensitivity method allowing the measurement of the non linear dielectric 
susceptibility of an insulating material at finite frequency. It has been developped for the study of 
dynamic heterogeneities in supercooled liquids using dielectric spectroscopy at frequencies 0.05 Hz< 
/ < 3 X 10^ Hz . It relies on the measurement of the third harmonics component of the current 
flowing out of a capacitor. We first show that standard laboratory electronics (amplifiers and 
voltage sources) nonlinearities lead to limits on the third harmonics measurements that preclude 
reaching the level needed by our physical goal, a ratio of the third harmonics to the fundamental 
signal about 10~^. We show that reaching such a sensitivity needs a method able to get rid of 
the nonlinear contributions both of the measuring device (lock-in amplifier) and of the excitation 
voltage source. A bridge using two sources fulfills only the first of these two requirements, but 
allows to measure the nonlinearities of the sources. Our final method is based on a bridge with 
two plane capacitors characterized by different dielectric layer thicknesses. It gets rid of the source 
and amplifier nonlinearities because in spite of a strong frequency dependence of the capacitors 
impedance, it is equilibrated at any frequency. We present the first measurements of the physical 
nonlinear response using our method. Two extensions of the method are suggested. 
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I. INTRODUCTION 



Measuring the nonlinear response of a physical system to an excitation is a way to in- 
vestigate physical properties often unreachable through the linear response. Understand- 
ing nonlinear effects allowed decisive breakthrough in condensed matter physics. Spin 
glasses 1^^, 3 
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Dinary 



3^ or 



40[ | effects in electrical transport, heating due to electric field 



excitation of supercooled liquids 4l|, |42| are a few among many topics where non linear 



measurements have proven to be a precious tool. 

The detection of harmonics of the fundamental response is a powerful method for studying 
nonlinear effects. It allows to get rid of the linear response signal, which usually is much 
larger than the sought nonlinear signals. We consider in this paper the nonlinear response 
of a dielectric system to a time-dependent electric field E{t). The method we present allows 
to extract very low level harmonics in the response to a sinusoidal excitation. The most 
general relationship relating the response (polarisation P{t)) to the excitation E{t) can be 
written as a series expansion in E (the even terms are forbidden because of the symmetry 
with respect to field reversal E{t) — * —E{t)): 

m 
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In this equation eo is the dielectric constant of vacuum, xi the hnear susceptibility and xs 

the cubic nonlinear susceptibility. The dots in Eq.[T] indicate an infinite sum involving higher 

order non linear susceptibilities xs? etc. Note that causality implies Xiit < 0) = 0. The 

Fourier transform of Eq. [T]for a purely a.c. field E = Eq cos{ujt) gives 
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where the polarization P and the susceptibilities Xi ^-^e now in the frequency domain and 



the dots indicate again infinite sums involving higher order terms. The response P{t) to 
E{t) = Eq cos(ti;t) can thus be written 

P(t)/eo = Re [{Eoxiiu;) + 3/4i?3x3(^) + -Oe-'"*] + Re [l/4Eo'x3(^)e-''"* + ...] + (3) 

where we have used the fact that because xi Xs ^ire real in the time domain, their 
Fourier transform verify Xii^) = Xi{~^) and X3(co'i, 0^2, cua) = X3{~^i, —^1) (the star 
denotes the complex conjugate), and the invariance of Xs by permutation of its arguments. 
For simplicity, we write Xsi^) = Xsi^y^y'-^) Xsi^) = Xsi'^^^'-^y'-^)- Eq. [3]can be written 

P{t) /eo = Eo{x'i cos ujt + Xi sin ujt) + S/AE^ix'^, cos oot 
+X5 sin cut) + ... + 1/ AEq {x's cos Scut + X3 sin Scot) + 

(4) 



where the susceptibilities Xi ^i-re written as a function of their real and imaginary parts Xi 
and Xi- For practical applications, the modulii and arguments \xi\ and 6i are rather used: 



P{t)/eo = ^0 Ixil cos(cut -6i) + 3/AE^ Ixsl cos^ - ^3) + 

+ ... + l/AE^ Ixsl cos(3c<;t - ^3) + •• 



(5) 



We see in the first and second terms of the rhs in Eqs [MS] that the nonlinear suscepti- 
bility X3 could be extracted from a measurement at the fundamental frequency by varying 
Eq. However, in most experiments ^ -^^o IXsl? thus the nonlinear part can hardly be 
separated from the much larger xi linear term. On the contrary measuring the harmonics 
yields directly the physical information contained in xs, X5; stc. In our case the physical 
information of interest is contained in xs (see section [Tll)- We thus chose to measure the 
third harmonics. However, its relative magnitude with respect to the fundamental was so 
low that we had to develop a special method for obtaining an accurate measurement. 

Experimentally, various methods have been used to extract a very small third, fifth, etc. 



harmonics signal 
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461]. A bridge technique has been often 



used: The first arm of the bridge contains the sample under study, and the second arm a 
well known impedance with zero nonlinear response j^]. This is the case for specific heat 
spectroscopy based on thermal diffusion into a thick sample from a thin metallic film that 



serves simultaneously as heater and thermometer 
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461] or similarly for the study of a 
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heating resistor where the third harmonics is related to its electrothermal parameters |35l. |45| . 
The bridge technique has also been used to measure the nonlinear inductive response of a 
superconducting film 25l| . the nonlinear behaviour of power piezoceramic materials 10| , etc. 
The bridge techniques evoked above cannot be used for our purpose because the balancing 
of the bridge at the fundamental frequency would be lost at the third harmonics: The strong 
frequency dependence of our sample impedance can hardly be mimicked by a combination 
of resistors and capacitors in the second arm. As a result, the third harmonics generated by 
the voltage source would not be cancelled in the bridge. 

In this paper, we report on a two samples bridge method allowing to measure properly the 
third harmonics (i.e. the third term in the rhs of Eqs[3}lS]), and eventually other harmonics of 
the polarization induced by the ac E-field excitation of a dielectric sample placed between 
the two electrodes of a plane capacitor. This is done through the measurement of the 
third harmonics of the current induced by the ac voltage excitation. Such a measurement 
yields directly the cubic nonlinear term xs of the response (assuming that the higher order 
contributions related to xs? X?) stc. to the third harmonics are weak which is the case in our 
experiment). We developped it for studying supercooled liquids such as glycerol. In part HTl 
we first briefiy summarize the physical interest of such measurements. Part IIIII is devoted 
to describing the experimental problems to solve. In part IIVI we present the two sample 
capacitors bridge method which allows to measure Xs- Part |V] gives our first measurements 
using this method. Finally, in part IVII we consider possible developments or extensions of 
this method. 



II. PHYSICAL MOTIVATION: DYNAMICAL CORRELATIONS IN GLASS 
FORMERS 



The physics of structural glasses still lacks a firm experimental basis for a growing 
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len the temperature T decreases towards the glass transition temperature 
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52| . A basic feature of glassforming supercooled liquids is the 



spectacular increase of the characteristic relaxation time Tq(T) as T decreases towards 
Tg. Tg is conventionally defined by Ta{Tg) ~ 100 s. is often obtained from dielectric 
50 . Ssl ]. Fig. [U shows an example of linear dielectric susceptibility Xi{^) 
measurement for glycerol (C3O3II8, Tg ~ 190 K) at T = 211.8 K. The imaginary part 



spectroscopy |48l. 
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Im{xi) is maximum, and the real part Re{xi) is approximately half its maximum plateau 
for Ua = iTxjTa. At present, the fast decrease of uJa with T has not received a unique 



microscopic physical interpretation, but a seminal concept [5^, |55| is that of cooperative 
effects: When Tg is approached, molecules belonging to larger and larger regions (called 
dynamical heterogeneities, DH) should move in a correlated way to allow relaxation. Corre- 
lation lengths probing such cooperative effects have been extracted ex per imentally, leading 

However 
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to length scales estimates of 4 — 12 molecular diameters 
these experiments were not able to test the expected increase of the correlation length as 
T decreases towards Tg which is of fundamental interest. For spin glasses the increase of 
the correlation length close to the critical temperature is associated to the divergence of Xi 
3)[l, y, y, ^, y, 0]. This suggests to investigate Xj(^) ^ 3) close to the structural 



z > 



glass transition, although no divergence is expected. Up to very recently, the non 
susceptibility Xsl"^) had never been measured in these systems (see however [41, 42 



mear 



6^) 



contrary to spin glasses. However, some important theoretical progress was made recently 
by Bouchaud and Biroli (GSJ who established that: 



60X1 ^j. njl 



(6) 



where is the molecular volume and 7i is a complex scaling function which should have a 
maximum (in modulus) for uot^ — 2tt. N*^^^ is the maximum of the number of correlated 
particles Ncorr{t) which increases with the time t, reaches its maximum N*^.^^. for t and 
then goes to for t — > 00. This is reflected in the frequency space by the fact that \l-L{ujTiy)\ 
should be maximum for uj uj^. Thus, measuring the nonlinear dielectric response Xs gives 
the size N*^^^ of the correlated regions and its T-dependence. 

What is the value of x^{'^) that can be expected from Eq. [6]? The answer cannot be very 
precise since little is known about Ti which, theoretically, should reach a maximum of "order 

63j. A conservative estimate of IxsLax = "^02;^ (1X3(^)1) for ^corr = 1 



1" for LOT "of order 1" 



can be drawn by assuming maXuj{\'H{ijj)\) = 1. We consider the case of glycerol at 200 K. 
Xi(0) should be replaced by Xi(0) — Xi(oo) = Ae ~ 72 {xi{^) = ^{^) ~ 1) because only the 
contribution of the molecular motion to the dielectric response is considered, ~ 0.115 nm^ 
is obtained from the density at 200 K [64'] and the molecular mass 92.09 g. As a reference, 
we take |e(t^„ = 2n/Ta)\ = |1 + Xii^a)\ ^ 52, and find by using Eq. [6], \x3\max I - 
3.685 X 10^^'' m^.V~^. We shall use this value in section UlI Al to determine the sensitivity 



required in our experiments. 



III. EXPERIMENTAL SETUP AND REQUIREMENTS FOR THE ELECTRON- 
ICS 

A. Cryogenics and nonlinear dielectric measurements 

The experiments were performed in a cryostat connected to a cryogenerator with a base 
temperature of 10 K. The experimental cell is a closed metallic box placed in vacuum, 



related to the cold stage of the cryogenerator through a thermal impedance [53|]. The 
temperature T in the cell is set by a PID LakeShore® 331 controller which regulates the 
heating power flowing through the thermal impedance. The cell contains two independent 
plane capacitors of equal surfaces S but different thicknesses Ltuck and Lthin (The reason 
for two capacitors instead of one is given in section II Vp . The capacitors are used to 
measure the dielectric susceptibilities of a glassforming liquid (glycerol) placed between the 
electrodes. For each capacitor, the two electrodes are immersed in the liquid. We used 
two experimental setups A and B. In setup A, the electrodes of the two capacitors are 
polished and gold plated brass squares (5* ~ 5.5 cm^), separated by three 0.03 cm^ Mylar® 
discs of thickness Lthin — 30 /xm for the thinner sample and Lfhick — 60 fim for the thickest 
sample. Setup B is an improved version, in which the electrodes are gold plated copper disks 
(metallic mirrors) with Lthin — 19 /im, Lthick — 41 fim and S = 3.14 cm^. We minimized 
the spacers volume in order to minimize their contibution to the dielectric response of the 
capacitors. The supercooled liquid was allowed to flow in and out of the volume between 
the electrodes in order to avoid pressure effects due to the different dilatation coefficients of 
the supercooled liquid and the spacers. The results presented in sections IIIH HVl (resp. |V|) 
were obtained using setup A (resp. B). A pressure of 2 bars of Argon was set in the cell at 
room temperature to ensure that, at the working temperature of ~ 200 K, the pressure in 
the cell remains above 1 bar, preventing the formation of bubbles of the gas adsorbed at the 
surface of the electrodes or dissolved in glycerol. Coaxial shielding was ensured all along the 
circuit, from the measurement apparatus at room temperature down to the experimental 
cell at low temperature. 

When a voltage V(t) = Vq cos(co't) is applied to the electrodes of a plane capacitor of 
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thickness L and surface S", the resulting field E{t) = V{t)/L induces a polarization P(t) 
and an electrical displacement D{t) = eoE{t) + P{t). As a result, a current I{t) fiows in 
the circuit, the current density I/S being the time derivative of D{t). Using Eq[3l in which 
only the first terms of the fundamental and third harmonics series are kept, we thus have 



where and I{3uj) are complex numbers giving the magnitude and the phase of the two 
components of the current. The linear part of the current is given by I{uj) = Y{uj)Vo where 
the admittance Y{uj) = G + iuC comprises the conductance G = eoUjIm{xi{^))S/ L and the 
capacitance C = eoRe{l + xi{^))S/ L. As depicted in the upper inset of Fig. [H we measure 
the current through the voltage drop Va{uj) across a resistor r = 1 kQ in series with 
the capacitor. All the voltage measurements presented in this paper were performed with a 
standard commercial phase sensitive lock-in amplifier with an input impedance Zl of 10 MQ 
in parallel with 25 pF. The accuracy of the Va{(^) measurement is typically not much better 
than 0.1%, which justifies that in the analysis of I{uj) we neglect the contributions of the 
nonlinear terms of the series giving the response at the fundamental frequency. The voltage 
sources were limited to Vs{uj) < 7 V (rms) for the measurements with setup A presented in 
sections llllt II VI and lIV Bl and to Vs{uj) < 14 V (rms) for those with setup B presented in 
section |Vl The frequency range was about 0.05 Hz< uj/2tt < 3x 10^ Hz. All the voltage 
and current magnitudes presented hereafter are rms. 

The calculation of the third harmonics component of the current, /(3ti;)e'^'^* (see Eq. [7]) 
from the time derivative of D (neglecting the contribution of the terms proportional to Xs, 
X7,--), gives 



In what follows, we investigate the measurement of I{3uj) (and I{uj)) for our liquid dielectric 
capacitors (that will be called "sample capacitors" or "samples" for simplicity). As the 
third harmonics current is simply added to the fundamental (Eq. [7]), the sample can be 
represented for what concerns the third harmonics by a current source given by Eq. [S] and 
placed in parallel with the liquid dielectric capacitor, whose complex admittance is Y{3uj) 
(see lower inset of Fig. [1]). Note that due to the Thevenin theorem, one could, as well, 
represent the sample at 3u! by a voltage source Vsampie = I{3u!) x Z{3uj) in series with the 
impedance Z{3uj) = 1/Y{3uj). 



lit) = Re{I{u)e-"^' + I{3u)e 



—i3Lxjt 



(7) 




(8) 
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From the estimate of Ixsl in section |TT1 we can determine the required sensitivity of our 
measurements. The complex currents I{uj) and /(Stu) being obtained from the complex 
fundamental and third harmonics components of the displacement, we have /(a;) = iujS{l + 
Xi)eo-^o and /(Scj) = Si/AuSxs^oEo, thus 



\I{uj)\ 4|l + xiM 



~ 2.76 X W'^^E^, (9) 



where the numerical factor is obtained 
according to the theoretical prediction [] 



br uj = Ua, and Eq is in V/m. We have assumed, 



63| that |x3(ti;a)| — Ixs I max' Vised the value of 



IXslmax estimated in section [TTl The maximum field Eq in our experiments is 200 kV/m for 
setup A (7 V on 30 fim), and 740 kV/m for setup B (14 V on 19 fim), thus \I{3lj)/ I{uj)\ 
~ 1.1 X 10^® or 1.5 X 10~^. The required sensitivity on the measurement of \I{3uj)/I{uj)\ is 
lower than these values, about 1 x 10~^ because of the uncertainty on Ti. and to allow the 
measurement of Xsi^) when its magnitude is below IxsLax- ^^^^ look for a setup whose 
relative parasitic contributions at 3a; remain below 10^^. 



B. Nonlinear behavior of the lock-in amplifier 

In this section we see why the simplest possible circuit, depicted in the upper inset of 
Fig. [H cannot work to detect I{3uj). For simplicity we assume that the resistance r is small 
with respect both to the impedance of the sample |2'(co')| and to the input impedance of the 
lock-in amplifier |Z/^(c<j)|. The voltage drop accross r at the fundamental frequency is thus 
Va{uj) ~ Vsr/Z{u). Fig. [2] shows that applying a lu signal directly at the input of the lock- in 
amplifier induces the measurement of a rather important 3uj signal. We found by using the 
method described in section IIII CI that this signal does not come from the source (which 
is that of the lock-in amplifier) which is expected to generate small harmonics in addition 
to the luj signal. It thus comes from the nonlinearities of the lock-in amplifier itself. We 
checked that the behavior reported in Fig. [2] does depend neither on the frequency, nor on 
the output impedance of the source tq. 

The magnitude of the 3uj signal due to the lock-in amplifier nonlinearities is much larger 
than the one expected from a glycerol sample in the simplest possible circuit, i.e. the sample 
in series with a measuring resistance r, see upper inset of Fig. [H Assuming r ^ |Z(c(j)|, one 
finds that /(3u;), see Eq. [HI yields a voltage VA,sampie{3uj) ^ rl{3uj), while I{uj) ~ Vs/Z{lj), 
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thus: 



(10) 



K Ziu) I{uj) 

We have seen in section UlI Al that a typical maximum expected value of |/(3co')//(a;)| 
was about 10~^ - 10~^. There is an optimum of the ratio r/\Z{uj)\ which maximizes 
\VA,sampiei^>^)\ (sec scctiou IIV A|) . and a reasonnable value is 0.1. Thus Eq. [10] gives 



VA,samp/e(3ci;)/K ~ 10^'' - 10^^. The corresponding Vs is of the order of 10 V and 
gives Va(co') ~ 1 V, thus by using Fig. [2] we find that the lockin nonlinearities give 
Va(3c(j)/Kj > 10^^, i.e. values much larger than the sample contribution. 

To summarize this section, we found that the nonlinearities of a standard lock-in amplifier 
are such that, due to our very low physical \I {"iuo) / 1 {uj)\ ratio, the luj component must be 
"removed from the signal", before it enters into the lockin for the Scj detection. This is 
why, in the next sections, we shall report results obtained with bridges performing a "lu; 
subtraction" . 

C. Nonlinear behavior of the voltage source 

1. Two sources bridge 

Let us consider the circuit depicted in the inset of Fig. [31 It uses two luj voltage sources 
Vs^i{uj) and Vs^2{'^) with a common ground and a tunable relative phase shift. The lock-in 
amplifier (input impedance Zl) is used to measure the voltage at point A which is connected 
to the reference resistor R and to the sample capacitor (impedance Z{uj)). Considering, for 
simplicity, the case R,\Z\ <^ \Zl\, we calculate the voltage at point A at the fundamental 



frequency : 



^^^"^^ = zuJ)Tr • ^^^^ 



The balancing condition which fulfills Va{uj) = (see section IIII Bp is 



(12) 



We used a Tektronix® AFG3102 dual channel voltage source for which the relative phase 
between K,i('^) and Kj,2(^^) is tunable with an accuracy of 0.01 degree, while |Kj i(co')| and 
1^5,2(1^)1 can be tuned with a relative precision of 10"''. As a consequence, |Va(co')/Vs,2(^^)| 
can be made as small as 5 x 10"^. As a result, we find by extrapolating the curve of Fig. [2l 



that the lock- in nonhnearities contribution to V^(3ct;) remains much below 2 nV for the 
maximum voltage K,2(ti^) = 7 V (which gives |VA(ci^)| ~ 5 x 10"^ x 7 V = 0.35 x 10"^ V). 

This value <^ 2 nV is small in comparison with the expected physical signal VA^sampiei'^^) ■ 
The latter can be estimated by using the circuit equivalent to the sample at Suj: a current 
source I{3uj) in parallel with the sample of impedance Z{3uj) (see section [Til Al and Fig. [T]). 



where we have assumed R, \Z\ ^ \Zl\- We chose R ~ 0.0841^(^7)1 (The maximum voltage 
cannot be applied to the sample if \R/Z{luj)\ > 1). For the capacitor we used {Lthm = 30 
/im), u = 2n/Ta (see section [TTl) and K,2(^^) = 7 V, and by using the value \I (Su) / 1 {uj)\ ~ 
10~^ obtained in section IIIl Al as the required sensitivity of our measurements, as well as 
\Z{ua,)/Z{3uj^)\ ^ 1.85, we find \VA,sampie{^^)\ ^ 6 x IQ-^ V > 2 nV . Thus the problem of 
the lock-in nonlinear contribution is solved by the two sources bridge. 

Despite this important feature, we shall see now that the two sources bridge does not 
reach the required resolution of |/(3co')//(ci;)| ^ 1 x 10^^. This comes from the fact that 
the two sources, as any active electronic device, have a non-zero harmonic distortion: They 
generate ^,1(3^^) and 1^,2(3^^) voltages in addition to K,i(^^) and 14, 2(^)- Thus a possible 
remaining problem is the contribution of Vs^ii^Su) and K, 2(3^1;) to V^(3u;). To reduce it, we 
low-pass filtered the sources outputs with a dual channel six poles elliptic active filter, with 
a corner frequency fc chosen close to the working frequency uj/27t. The filter damps the Slj 
component of the incoming signal, while it does not affect the lu component. We verified 
that this filtering reduced significantly the harmonics magnitude. Fig. [3] shows the measured 
1^(3^7) for the two sources bridge at u = Ua- The results are the same for the thin and 
the thick sample of setup A {Lthin = 30 /xm and Lthick = 60 yum), which shows that the 
measured signals do not come from the samples. Indeed, Eq. [H] gives a current |/(3c(j)| eight 
times larger for the thin sample than for the thick one. We took R ~ 0.084 |2'(u;)| in order 
to keep similar values of Vs^i{u)) and Kj,2(^) for the thin and the thick sample measurements, 
thus, from Eq. [131 the expected ratio of V^(3c(j) for the two samples is 4, instead of the "no 
change" result of Fig. [3l Thus we conclude that the third harmonics generated by the voltage 
sources are responsible for the results of Fig. [3l We verified that the expected physical value 
VA,samp/e (3ci;) (sec scctiou [Til Ap is indeed lower than the measured Va{Suj). 

A closer look at the equilibrium condition (Eq. [12]) of our two sources bridge reveals 




R/Z{uj) I{3uj) 
1 + R/Z{3uj) I{uj) 



(13) 
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that once the circuit is balanced at loo, it should not be in general balanced at 3uj, for two 
reasons. First, because of the strong change of Z{uj) when going from luj to (due to the 
Xi{^) dependence, see Fig. [T]) while R remains constant. Note that to mimick the Z{uj) 
dependence by using an impedance made of resistors and capacitors instead of R would 
be hardly feasible because its components should be changed and tuned for each frequency 
and temperature. Second, and more importantly, because K,i(3c<j) differs from Kj,2(3co') as 
can be deduced from the comparison of the signals measured when the two sources are 
exchanged (see Fig. [3]). 



2. Measurement of the harmonics signals generated by the sources 

The two sources bridge cannot be used to measure the nonlinear response of the samples, 
however it can be used to measure with a great accuracy the harmonics Vs^i{nuj) and ^,2(^1^) 
of our voltage sources. This is of general interest for characterizing the harmonic distortion 
of any source, and it will be used in section IIVI (Note that we propose another method 
in section |V]). We now consider only the thickest sample of setup A for which I{3uj) is 
neghgible. Assuming that the voltage measured at point A at frequency nu,VA{nu), is due 
to the sources, we find for the circuit in the inset of Fig. [3] 

VAinUj) = — r — r , (14) 

^ ^ y + Y{nuj) + YL{nuj) ^ ' 

with y = 1/R, Y{nuj) = 1/Z{nuj) and YL{nuj) = l/ZL{nuj). For any given n, another 
equation is necessary to determine the unknown complex quantities Vs^i{nuj) and Vs^2{nuj). 
It is obtained by exchanging the two sources in the circuit. The measured signal becomes 

VA,ex.- 

y + Y{nu!) + YilnLJ) 

Note that Vs^i{nuj) and K,2(^'^) are two functions of respectively ^^,1(1^^) and ^,2(1^^) 
which can be different from each other. Thus, for Eqs. [T4fT5] to contain the same unknown 
quantities, we had to choose R = \Z{lu!) \ in order that |K,i(l^^)| = |K,2(1'^)| because of the 
equilibration condition (see Eq. [T2]) . As a result Eqs [T^ and [T^ can be easily solved. 

Figure m displays the harmonics Vs^ii^nuj) of source 1 and K,2(^i^) of source 2 obtained 
by using the two sources bridge method, for uj/{2'k) = 88 Hz. The main features are: i) The 
second harmonics of both sources are similar in magnitude and phase. |V^_j(2a;)| increases 
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as the power two of the lu voltage. ii) The third harmonics is twice smaller for source 2 
than for source 1. One finds 1^4,2(3^7)1 oc |\4,2(1'^)P'^, i-e. a nearly cubic dependence which 
could have been mistaken with the physical signal if the tests mentionned above had not 
been performed. Hi) The Total Harmonic Distortion (THD), given by summing up all the 
harmonics magnitudes n > 2 and dividing by the magnitude of the fundamental signal, is 
dominated by the second harmonics. Thus, from Fig. HI at K = 7 V the THD is of the 
order of 5 x 10^^, well below the 10~^ specification which is usually ensured by high quality 
electronic devices. These features remain basically true for all the frequencies we studied 
in the 1-100 Hz range. Fig. Ogives the second and third harmonics of the second source at 
4.28 Hz that will be used in the next section. 

We summarize this section IIIII by emphasizing the two main requirements that have to 
be met to measure the nonlinear response of our samples: First, the luj part of the signal 
has to be suppressed before amplifying the sought 3u! signal, and this can be done by using 
a bridge technique. Second the bridge must be balanced at u and 3uj. 

IV. A SETUP ALLOWING THE MEASUREMENT OF X3(w) OF SUPERCOOLED 
GLYCEROL 

A. One source bridge with two samples 

Let us consider the circuit depicted in the inset of Fig. O It uses only one source and 
is inspired from the Wheastone bridge, but it contains two liquid dielectric capacitors of 
different thicknesses Lthm o-nd Lthick- The right arm of the bridge contains the thin capacitor 
("thin sample") of impedance Zthin oc Lthin in series with a chosen impedance Zthin, while 
the left arm contains the thick capacitor ("thick sample") of impedance Ztuck oc Ltuck in 
series with a chosen impedance zthick- The signal Vm is measured between the two middle 
points of the two arms with our lock-in amplifier in differential mode. The input impedances 
of the lock-in are much larger than \zthin\ and \zthick\, thus the measured voltage Vm,s at any 
frequency, due to the source voltage Vg at the same frequency is 



where the impedances are taken at the frequency considered. The bridge is balanced for 




(16) 



^thin^ thick ^thick^thin- 



(17) 
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The key point is that if we choose Zthin and Zthick of the same nature (e.g. two resistances, or 
two capacitances), once Eq. [T7]is fulfilled at a given frequency, it is fullfilled at any frequency. 
This is because the relative variations of the two samples impedances with frequency are the 
same on both sides of Eq. [T71 As a result, if the equilibration condition (Eq. [T7|) is fulfilled 
at luj, all the harmonics Vs{nuj) generated by the source will give a zero contribution to Vm- 
This two-capacitors bridge thus allows to get rid of both the nonlinearities of the amplifier 
(see section UlI Bp and of the harmonics generated by the source (see section UlI C|) . 

Obviously, the physical contributions from the samples yield a non-zero measured signal 
Vm,sampie{3^) bccausc the two currcut sources at 3u! corresponding to the samples verify 
Ithini^^^) oc L~i^^^ and Imcki^uj) oc L~^^^^ (see section [ITS])- For a typical Lthick ~ '^Lthin, 
the factor 8 between the two currents is not compensated by the factor 2 between Zt^ick and 
Zthin'- The relation between the measured signal and the two physical currents is 



where Yttin = I /Zthin, Ymck = I /Zthick, Vthin = l/zthin and y thick = l/zthick- By using 
EqsEHZIand Ythick/ymn = L thin / L thick, Eq. [H becomes 



In practice, zthin and zthick were resistors {rthin and r thick), but we improved the equilibration 
by adding a small capacitor Cthin in parallel with rtun- Cthin compensates for the unavoidable 
stray capacitances between the circuit and the ground. The equilibration was realized by 
tuning rthin, r thick and Cthin to minimize |Kn| at loj. In practice, the equilibration at noj {n 
> 1) is not perfect: Small asymetries between the two samples, stray capacitances which 
are not proportionnal to the samples thicknesses etc make that Vm{nuj) is not exactly zero. 
The bridge is thus characterized by the quality factors p{nuj) = \Vm{nLj)/Vs{nuj)\, which 
can be measured by tuning the source fundamental frequency at nuj. With setup A (S* ~ 5.5 
cm^, Lthin — 30 /im, Lthick — 60 /im, see section ITII Ap . we reached e.g. p{luj) ~ 3 x 10^^, 



p{2uj) ~ 5.9 X 10-^ and p{3uj) ^ 1.2 x lO^^ ^/27r = 27 Hz and u; ^ (T = 208 K). With 
setup B (S ~ 3.14 cm^, Lthin — 19 pm, Lthick — 41 pm), we obtained p{luj) ^ 2.7 x 10~^, 
p{2uj) ~ 6.6 X 10-^ and p{3uj) ~ 1.4 x 10"^ for uj/2tt = 29 Hz and ~ cj„ (T = 209.5 
K). Note that smaller values of p{3uj) are obtained in the range u <^ uja (i.e. when the 




(18) 



m,sample\3U! 



ythin{3uj) + Ythin{3uj) ythick{3uj) + Ythicki^^) ' 




(19) 
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capacitances of the samples are the largest): for example, with setup B, for T = 204. S/C and 
u/ua = 0.0567, we find p{3uj) ^ 2 x 10~^ while p{luj) is only one order of magnitude lower. 

There is an optimum of the ratio rthin/ \Zthin{^^^)\ = r thick / \ZthickO-^)\ which max- 
imizes the measured voltage \Vm,sampie{'i^)\- If this ratio decreases, the voltage Vq = 
Vsiluj) \Zthinl ijthin + Zthin)\ applied to the samples increases, thus |/(3co')| oc increases 
too; but at the same time, the contribution of I{?>uj) to Vm,sampie{^>^) decreases because 
it results from this current flowing through rthin (or rtuck) in parallel with Zthini'i^) (or 
Zthicki'i^))- We found that this optimum was rthin/ \Zthin{^^)\ = rthick/ \Zthick{^^)\ - 0.36. 
This value depends only weakly on the frequency and on Ta{T) (see section HTl) . 



B. Detection of the nonlinear dielectric response of glycerol 

We show in this section that by using the two sample bridge with setup A (see sec- 
tion [TllA]) the source 2 (whose harmonic distortion has been studied in section UlI Cp . 



we reached the sensitivity needed to detect the physical signal. We compare the measured 
voltage due to the source harmonics Vm,s{'i^) to the expected physical signal Kra,samp«e(3co'), 
for the maximum source voltage K,2(t^) = 7 V at uj /(2v:) = 4.28 Hz (tu/a;^ ^ 1 for T = 203.7 
K), with setup A {Lthin = 30 fim, Ltuck = 60 /im). Using the results reported in Fig.[5]and the 
measured p(3cj) ^ 1.4 x 10^^ we find |V;„,,(3cj)/V;,2(MI ^ 7 x 10"^ By using \I{3uj) / I{uj)\ 
^ 1.1 X 10^^ (see section UTTAl and Eq|9]) and Eq. [19], we obtain |Ka,samp«e(3co')/ys^2(lt^)| ~ 
8 X 10^*^. The measured voltage is Vmi/i^jj) = Kn,s(3Lij) + Vm,sampiei^>^) i and the physical 
signal should contribute to Vmi/ioj). 

Figs [6] and [7] give the measured Vmi^uj) and Kn(3c<j) for the same frequency 4.28 Hz and 
temperature 203.7 K. Fig. [6]shows that the measured V^(2u;) can be calculated by assuming 
that it is only due to Vs{2uj) generated by the source, in agreement with I{2uj) = due to 
the field reversal symmetry mentionned in the introduction. On the contrary. Fig. [7] shows a 
discrepancy between ^(3^^) and the Vm,s{'^^) value that should be measured if it was only 
due to ^(3^1;) generated by the source, strongly suggesting that a physical signal is present. 
As for n = 2,we used 1^^,8(3^7)1 = p{3uj)\Vs{3uj)\, with a proper treatment of the phase. The 
physical signal should grow as Vs{uj)^, while Vs^2{3uj) grows as Vsi^uY'^ (see section 1111 C 21) . 
This difference of exponents explains that Vm{3uj) grows slower at high voltages where the 
physical signal should dominate progressively Vm,s{3uj). This also explains the evolution of 
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the phase. The Vm{3u!) vs. Vs{uj) data can be fitted by a sum of ^71,5(30;) and a physical 
signal A[K(i^)]^ where A is a complex number. The same can be done at other temperatures 
and frequencies. However, we rather present in the next section the results obtained with 
setup B using another source: The improved quality factor p and |/(3ti;)//(co')| ratio for this 
setup, and the lower THD of the source allowed the first accurate measurement of Xs- 



V. MEASUREMENT OF THE NONLINEAR DIELECTRIC RESPONSE OF A 
SUPERCOOLED LIQUID 

We present in this section our first measurements of the nonlinear dielectric susceptibility 
of supercooled glycerol. We used the two samples bridge, with setup B (see section IIII A[ 
Lthin — 19 /im, Lfhick — 41 /um and S" = 3.14 cm^). We used a voltage source with a lower 
total harmonic distortion (THD) than the sources used in sections IIIH IIVI It was a SRS® 
DS360 ultra low distorsion function generator, with < 14 V and a typical THD ranging 
from -98 dB to -109 dB (for Vs = 10 V) in our frequency range. Note that due to the very 
low level of our physical signal, our bridge technique remains necessary even with such a 
low THD. Fig. [8] shows the current /t/ij„(3c<j) calculated from the measured Kn,samp/e (3lij) by 
using Eq. [191 as a function of the source voltage magnitude \Vs\, for T = 210.2 K, uj/2'71 = 
uja/27r = 43.76 Hz. A clear \Ithin{3uj)\ oc [^^(to')!^ dependence is found, while the phase of 
I{3uj) does not depend on |\^5(co')|. These results are a strong indication that the origin of 
the signal is really the nonlinear response of the samples. A constant phase suggests that no 
other signal than the physical one is present, contrary to the results of the previous section. 
The magnitude is of the order of the estimated physical value |/(3ci;)//(co')| ~ 1.5 x 10~^ for 
Vs = UV (see section [ml) which gives \I{Suj)\ ^ 3 x 10"^° A. 

To confirm the physical origin of our measurement, we compared these results to the data 
obtained by using a passive notch filter method. The principle of this method is depicted 
in the inset of Fig. [8j The source voltage Vs{uj) is apphed to a sample capacitor (impedance 
Z) in series with a resistor r. The voltage a point A (Va(3u;)) is measured by the lock- 
in amplifier. However, contrary to the method depicted in section IIII Bl the fundamental 
component of the voltage at point A, Va{uj) is strongly attenuated by a "twin-T" passive 
notch filter with a center frequency fc = u;/27r. The twin-T notch filter is made of three 
capacitors C, C, 2C and three resistors 2R, 2R R, and fc = 1/2ttRC 65|. A passive filter 
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avoids the nonlinearities related to active components, and we verified that the resistors and 
capacitors did not yield noticeable 3uj contribution to Va- The attenuation coefficient at lu 
(a few 10"^) was such that the lu component at the lock-in input yielded a negligible 3lo 
harmonics (see Fig. [2]). To subtract the contribution to Va(3u;) of the third harmonics due 
to the source Vs{3uj), the measurement was repeated with a circuit in which the sample was 
replaced by a resistor of impedance |Z(ti;)|. Finally, the currents I{3uj) obtained by using 
this method with the thin or the thick capacitor were proportional to the thicknesses at the 
power 3 as expected (see Eq. [8]). The current I{3uj) obtained with the twin-T ffiter method 
is in remarkable agreement with the current obtained using the two sample bridge, both in 
what concerns the phase and the magnitude (see Fig. [H]). 

We also confirmed the purely physical origin of our results of Fig. [S] by verifying that the 
two possible contributions to Vm{3u) due to the source and to the lock- in amplifier were 
negligible in the case of the two samples bridge. As explained above, Vs{3uj) was measured 
by using the twin-T notch ffiter method with the sample replaced by a resistor. Its value 
depends on uj and Vg, but we found that its contribution to Vm{3uj) remained much smaller 
than the physical contribution for the two samples bridge, e.g. for T = 210.2 K and uj/{2t[) 
= 43.76 Hz, it is 3 orders of magnitude below. The contribution of the lockin was negligible 
because the equilibration of the bridge led to Vm{uj) values below 1 mV: Fig. [2] shows that 
the resulting third harmonics is well below 40 nV, a negligible value with respect to the 
measured Vm.{3uj) of a few fiV. This remained true when the lockin-in amplifier was used in 
the differential mode provided that the voltages on each input remained below about 1.4 V. 
Another possible spurious contribution to the measured third harmonics could come from 
the resistors rthin, f thick and the capacitor Cthin- We verified in experiments using bridges 
made only of resistors or of capacitors or of combinations of both that the nonlinearities of 
metal ffim resistors and polystyrene film capacitors (we rejected ceramic capacitors) were 
negligible. Finally, another verification of the physical origin of our data was that |Kn(2co')| 
remained much weaker than |Kn(3ti;)|. 

VI. POSSIBLE DEVELOPMENTS OF THE TWO SAMPLES BRIDGE METHOD 

The equilibration (Eq. [T7]) of the two-sample bridge is realized by tuning Ztun and Ztuck 
(two resistors in our case) in series with the two impedances of the samples, Ztun and Ztuck 
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(see section IIVI and the insets of Figs. [6|7|) . This is time consuming because in order to 
keep the ratio rthin/ \Zthin{,^^)\ = fthick/ \Zthick{,^^)\ at its optimum value, it has to be done 
for each temperature or frequency change. A possible solution is a four samples bridge: 
In such a device, the two sample capacitors of different thicknesses are kept, but the two 
impedances Zthin and -Zthjcfc are replaced by two other sample capacitors. Thus, neither the 
equilibration, nor the ratio of impedances in each arm depend on T or uo. The first (resp. 
second) arm of a four sample bridge would contain sample 1 (resp. 2) connected to the 
source and sample 3 (resp. 4) connected to the ground. The thicknesses are Lj {i = 1,4) 
and the electrode surfaces Si {i = 1,4). If the surfaces are the same in an arm (5*1 — 5*3 or 
5*2 = ^4), the contribution of the two physical currents (11(30;) and l3{3uj), or l2{Suj) and 
I^i^Suj)) to the measured 3u voltage Vm^sampiei'^^^) is zero. This is because the electric field at 
luj in the two capacitors in an arm would be the same: The different thicknesses (see Eq. [8]) 
are compensated by the voltage division between the two capacitors. Thus Si 7^ 5*3 and S2 
7^ 5*4 are needed. The balancing condition Z1Z4 = Z2Z3 leads to L1L4S2S3 = L2L3S1S4. 
Following the lines of section IIV Al the measured Sou voltage is given by 

Ki,samp;e(3(^) S1S3L1L3 ^ q2\ S2S4L2L4 ^ q2\ frfr,\ 

K3(M (LiSs + LsS.r^^' {L2S4 + L4S2Y^^' ^ 

The two (or four) sample capacitors bridge methods can be generalized to any physi- 
cal situation in which a sinusoidal excitation F{t) = /2e[-Foe*'^*] gives a response R{t) = 
i?e[xi(a;)Foe*'^* + |x3('^)-^o^c^*'^* + •••] where xs, Xs are to be measured. If Fq depends on a 
parameter A analogous to our capacitor thickness, it is possible to design a bridge in which 
each arm contains an experimental unit with two different values of A. Again, the advan- 
tage of such a bridge with respect to a bridge where the sample response is balanced by the 



response of an impedance 
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45| (or possibly with two sources as discussed 



in section IIII Cp is that if Xi depends on the frequency, the equilibration at lu implies also 
an equilibration at nu. Among the possible applications of this principle, we can think of a 
bridge devoted to the measurement of the non linear magnetic susceptibility: The capacitors 
of our bridge would be replaced by coils with sample cores, and the A parameter would be 
the number of turns of the coils. 
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VII. CONCLUSION 



Many experiments in physics consit in applying an "excitation" (electric, magnetic, me- 
chanical, etc.) to a sample at a finite frequency and measuring the response. The nonlinear 
part of this response is present in its harmonics. Their measurement is a hard task for 
weak nonlinear responses because the electronic devices used for the excitation and for the 
measurement are always non linear at a certain level. Our experiment, devoted to the 
measurement of the nonlinear response of dielectric supercooled liquids, measures the third 
harmonics component of the polarization when an electric field is applied to a sample layer 
in a plane capacitor. It consists in measuring a current /(Scu) at a frequency three (pos- 
sibly five, etc.) times that of the linear response I{uj) when a voltage Vs{uj) is applied to 
the capacitor. Both the lock-in amplifier and the voltage source yield a third harmonics 
component which precludes measuring the nonlinear response in a simple experiment such 
as this presented in section IIII B[ To overcome this problem, we have used a method based 
on a bridge with two capacitors of different thicknesses. We have shown that such a device 
strongly reduces the two spurious third harmonics component of the measured signal: That 
of the voltage source and that of the lock-in amplifier. This "two-sample" bridge, used with 
a low distortion voltage source allowed to reach our goal of a sensitivity better than 10~^ 
and to realize our first measurements of the nonlinear response of a supercooled liquid close 
to the glass transition. 
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FIG. 1: Measurement of the dielectric susceptibility e = e' + ie" = xi + 1- The experiment is 
shown in the upper inset: The capacitor with glycerol as dielectric material is equivalent to a 
non-dissipative capacitor C{uj) oc e' in parallel with a conductance G{uj) oc ue". The impedance 
Z{uj) = {G + iCLo)~^ is obtained by measuring the voltage drop across the resistor r at point A 
when the excitation voltage V = FsCoswt is applied in B. The measured values of C and G/lo 
are shown as a function of uj/2tt for T = 211.8 K. In the lower right corner, the equivalent circuit 
corresponding to the physical third harmonics source is shown: an ideal current source in parallel 
with a capacitor C{3lo) and a conductance G{3u). 
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FIG. 2: Magnitude of the measured third harmonics voltage as a function of the source voltage 
when the latter is directly applied at the input of the lock- in amplifier (for a;/(27r) = 88 Hz). 
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FIG. 3: The third harmonics voltage |Va(3u;)| measured at point A vs. the source voltage |Kj(a;)| = 
max{\Vs^i{lu!)\ , |ys,2(l'^)|) applied to the sample for the two sources bridge depicted in the inset. 
The components of the circuit are the sample (impedance Z), the two voltage sources (V^^i(cj) and 
Vs^2i^))i the resistor R and the lock- in input impedance Z^. Closed diamonds with short dashed 
line: Source 2 applied to a thin sample {Lt^in = 30 fim). Open squares with continuous line: Same 
circuit, but with a thick sample {Lfhick = 60 ^um). Open triangles with long dashed line: Source 1 
is applied to the thin sample, and source 2 to R. The temperature of the samples is 211.8 K and 
a;/(27r) = u;a/(27r) = 88 Hz 
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FIG. 4: Magnitude (left axis) and phase (right axis) of the second and third harmonics generated 
by the source 1 (closed symbols) and source 2 (open symbols) of the two sources bridge, obtained 
by solving Eqs. ll4|15l for a frequency of 88 Hz. 
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FIG. 5: Magnitude (left axis) and phase (right axis) of the second and third harmonics generated 
by the source 2 of the two sources bridge, obtained by solving Eqs. I14|15l for a frequency of 4.28 
Hz. 
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FIG. 6: Magnitude (left axis) and phase (right axis) of the second harmonics Vra{2uj) measured in 
the two samples bridge depicted in the inset. The experimental data (closed symbols) are com- 
pared to the values (open symbols) calculated by assuming that 14n(2u;) is only due to the second 
harmonics generated by the source Vs{2uj) which was obtained from the two sources measurements: 
\Vm{2uj)\ = p{2uj) \Vs{2uj)\. uj/{2tt) = 4.28 Hz and T = 203.7 K 
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FIG. 7: Magnitude (left axis) and phase (right axis) of the third harmonics 14n(3u;) measured in the 
two samples bridge depicted in the inset. The experimental data (closed symbols) are compared 
to the values (open symbols) calculated by assuming that Vmi'i^^) would be only due to the third 
harmonics generated by the source V's(3a;). lo/{2tt) = 4.28 Hz and T = 203.7 K 




FIG. 8: Open circles: The current |It/ij„(3cj)| calculated from the measured Vm,sampie{'^^) in the 
two samples bridge, as a function of the source voltage magnitude |ys(w)|, for T = 210.2 K, uj /2Tr 
= 43.76 Hz. Closed diamonds: [/fft,m(3a;)| obtained using twin-T filter method. The continuous 
line is the power 3 dependence. Inset: The twin-T notch filter method (see text). 
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